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ABSTRACT 

Recent studies suggest that pulsars could be strong sources of TeV muon neutrinos 
provided positive ions are accelerated by pulsar polar caps to PeV energies. In such 
a situation muon neutrinos are produced through the delta resonance in interactions 
of pulsar accelerated ions with its thermal radiation field. High energy gamma rays 
also should be produced simultaneously in pulsar environment as both charged and 
neutral pions are generated in the interactions of energetic hadrons with the ambient 
photon fields. Here we estimate TeV gamma ray flux at Earth from few nearby young 
pulsars. When compared with the observations we find that proper consideration of 
the effect of polar cap geometry in flux calculation is important. Incorporating such 
an effect we obtain the (revised) event rates at Earth due to few potential nearby 
pulsars. The results suggest that pulsars are unlikely to be detected by the upcoming 
neutrino telescopes. We also estimate TeV gamma ray and neutrino fluxes from pulsar 
nebulae for the adopted model of particle acceleration. 
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1 INTRODUCTION 

Probable candidates of high energy neutrino radiation in- 
clude Gamma ray bursts, Active galactic nuclei etc (Halzen 
& Hooper 2002). Recently Link and Burgio (in this arti- 
cle we hereafter refer them as LB) (Link & Burgio 2005, 
2006) have shown that pulsars also could be strong source 
of high energy neutrinos. As per their results pulsars with a 
magnetic moment component antiparallel to the spin axis, 
which is expected in half of total neutron stars, could emit 
TeV muon neutrinos with fluxes observable by the operating 
or planed large area neutrino observatories. As a conjecture 
to be verified by observations LB considered that protons or 
heavier ions are accelerated near the surface of the pulsar 
by the polar caps to PeV energies. When these accelerated 
ions interact with the thermal radiation field of pulsar the 
A resonance state may occur provided their energies exceed 
the threshold energy for the process. Though radiation losses 
limit the maximum energy that can be attained by a nuclei 
in the acceleration process but such an energy condition is 
expected to satisfy for several pulsars. Muon neutrinos are 
subsequently produced from the decay of A particles. 

The presence of a hadronic component in the flux of 
pulsar accelerated particles should result in emission of high 
energy neutrinos and gamma rays simultaneously as both 
charged and neutral pions are produced in the interactions of 
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energetic hadrons with the ambient photon fields surround- 
ings the acceleration region. Here we would show that for 
the model adopted by LB, the estimated TeV gamma ray 
fluxes from several nearby pulsars are higher than the ob- 
served upper limits of fluxes. In the quest for reasons of such 
inconsistency we note that an implicit assumption in the LB 
estimation of neutrino flux is that the polar cap area is equal 
to the neutron star surface area. Such an assumption seems 
unreasonable in view of the polar cap geometry (Beskin, 
Gurevich & Istomin 1993). When this fact is taken into con- 
sideration, the stated inconsistency between the estimated 
and the observed gamma ray fluxes is found to disappear. 
In view of this observation a revised estimates of the neu- 
trino fluxes from few known gamma ray pulsars are obtained 
incorporating polar cap geometry. 

A young neutron star is generally encircled by pulsar 
wind nebula. Positive ions, after gaining energy from polar 
gaps will move away from the pulsar practically along the 
open field lines and will finally inject into the nebula. It 
is very likely that these energetic ions would be trapped 
by the magnetic field of the nebula for a long period and 
consequently they should produce high energy gamma rays 
and neutrinos by interacting with the matter of the nebula. 
We, therefore, estimate the expected flux of TeV gamma rays 
from pulsar nebulae and by comparing with the observation 
for a couple of well known nebulae we check the consistency 
of the model. We also calculate the flux of TeV neutrinos 
from couple of pulsar nebulae. 

The organization of the article is as follows. In the next 
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section after describing generation of the TeV gamma rays 
and neutrinos in a pulsar environment, the expected flux of 
gamma rays from a pulsar is estimated for the polar cap 
model of acceleration as used by LB. By comparing the 
model predicted 7 ray fluxes from some potential young pul- 
sars with the observations, importance of inclusion of polar 
cap area in the calculation has been stressed in section 3. In- 
corporating such a feature (revised) event rates in a neutrino 
telescope at Earth due to few potential pulsars are obtained 
in section 4. In section 5, fluxes of TeV gamma rays and 
neutrinos from pulsar nebulae are obtained for the adopted 
model and finally the results are concluded in section 6. 



2 TEV GAMMA RAYS AND NEUTRINOS 
FROM PULSARS 

Several detailed mechanisms have so far been suggested for 
acceleration of particles by pulsars those include the popular 
polar gap (Ruderman $ Sutherland 1975; Arons & Scharle- 
mann 1979; Daugherty & Harding 1996; Harding & Mus- 
limov 1998) and the outer gap models (Cheng, Ho & Rud- 
erman, 1986). In the former model, acceleration of particles 
takes place in the open field line region above the magnetic 
pole of the neutron star whereas in the case of outer gap 
model it occurs in the vacuum gaps between the neutral 
line and the last open line in the magnetosphere. Thus, the 
region of acceleration in the polar gap model is close to the 
pulsar surface, while the same in the outer gap model is close 
to the light cylinder. 

In the polar cap acceleration model, particles are ex- 
tracted from the polar cap and accelerated by large rotation- 
induced electric fields, forming the primary beam. The max- 
imum potential drop that may be induced across the mag- 
netic field lines between the magnetic pole and the last field 
lines that opens to infinity is A(f> = BsR 3 Q, 2 /2c 2 (Goldreich 
& Julian, 1969) where Bs is the strength of magnetic field 
at neutron star surface, R is the radius of the neutron star, 
Q is the angular velocity and c is the speed of light. Accord- 
ingly, for young millisecond pulsar with high magnetic fields 
(Bs ~ 10 12 Gauss), the magnitude of the potential drop 



could be as large as 7 x W 1S B 12 P~; volts (B s = B 12 x 10 12 
Gauss). But in young pulsars, the electric field along mag- 
netic field lines is likely to be screened well below the vacuum 
potentials by the onset of electron-positron pair cascades in 
the strong magnetic field (Cheng & Ruderman 1977) and 
this would limit the potential to about 10 12 eV. Though the 
flux of synchrotron radiation observed from the Crab and 
other pulsar wind nebula (PWN) indiactes such a possibil- 
ity but questions like where in the magnetospher are pairs 
created and how many are created, are still not settled and 
one cannot rule out the possibility that ions can reach en- 
ergies equal to a significant part of the total potential drop 
through polar cap acceleration, particularly in view of the 
observation of gamma radiation of energies above tens and 
even hundreds of TeV from Crab and other PWN. 

LB (Link & Burgio 2005, 2006) conjectured that pro- 
tons or heavier ions are accelerated near the surface of a 
pulsar by the polar caps to PeV energies (corresponds to 
no/little screening) when /i.fi < (such a condition is ex- 
pected to hold for half of the total pulsars). When pulsar 
accelerated ions interact with the thermal radiation field 



of pulsar the A resonance state may occur provided their 
energies exceed the threshold energy for the process. The 
threshold condition for production of A resonance state in 
P7 interaction is given by 



e p e 7 (l — cosOp-f) > 0.3 GeV 



(1) 



where e p and e 7 are the proton and photon energies re- 
spectively and dp-y is the incident angle between the proton 
and photon in the laboratory frame. In a young neutron star 
with surface temperature T x the energy of a thermal pho- 
ton near the surface of the neutron star is 2.8fcT oc ,(l + z g ), 
z g (~ 0.4) being the gravitational redshift. This implies that 
in a young pulsar atmosphere the condition for production 
of the A resonance is Bi2-Pm 2 To.i kev 1? 3 x 10~ 4 (Link 
& Burgio 2005, 2006) where T .i keV = (kT^/O.lkeV), 
Too ~ O.lkeV being the typical surface temperature of young 
pulsars. Such a condition holds for many young pulsars 
and thus A resonance should be reached in pulsar atmo- 
sphere. Gamma rays and neutrinos are subsequently pro- 
duced through the following channels: 



P + l 



P + 7T° -> P + 2 7 

nir + — > n + e + + v e 



Since the charge-changing reaction takes place just one- 
third of the time, on the average four high-energy gamma- 
rays are produced for every three high-energy neutrinos (or 
for every two high energy muon neutrino and antineutrino) 
when a large number of such reactions occur. 

The flux of gamma rays and muon neutrinos from pul- 
sars can be estimated as follows: The charge density of 
ions near the pulsar surface is p q ~ eZn , where n (r) = 
B 3 7? 3 r2/(47rZecr 3 ) is the Goldreich- Julian density (Goldre- 
ich & Julian, 1969) of ions at radial distance r. For accel- 
eration to take place there must be a charge depleted gap 
(here polar gap) and the densiy in the gap may be written 
as fd(l — fd)n , where fd (< 1) is the depletion factor which 
is a model dependent quantity (f d = corresponds to zero 
depletion). The flux of protons accelerated by a polar cap is 
therefore 



c/d(l - fd)n A pc 



(2) 



where A pc — rjAirR 2 is the area of the polar cap, r) is 
tha ratio of polar cap area to the neutron star surface area, 
which is taken as unity by LB in their work (i.e. in original 
calculation by LB A pc was taken as equal to the surface area 
of the whole neutron star). The canonical polar cap radius 
is given by r pc = R(nR/c) 1/2 (Beskin et al. 1993) and thus 



r, = OR/ (4c) 



(3) 



The protons accelerated by polar cap will interact with 
the thermal radiation field of the neutron star. The tem- 
perature of polar caps is expectedly higher than the surface 
temperature of neutron star but the contribution of polar 
caps on the thermal radiation field of a neutron star should 
be negligible because of their small surface area in compar- 
ison with the surface area of the whole neutron star. For a 
young neutron star with surface temperature Too the pho- 
ton density close to the neutron star surface is n 7 (i?) = 
(a/2.8fc) [(1 + z 9 )T 0O ] 3 , a being the Stefan-Boltzmann con- 
stant. Numerically n 7 (i?) ~ 9 x W 19 Tq A keV . At radial dis- 
tance r, the photon density will be n 7 (r) = n 1 (R)(R/r) 2 . 
The probability that a PeV energy proton starting from 
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the pulsar surface will produce A particle by interacting 
with thermal field is given by (Link & Burgio 2005) P c = 

1 - Ir P ^ where dP l P = -n~i{r)<?vidr . The threshold 
energy for production of A resonance state in p-y interaction 
as given by Eq.(l) increases rapidly with distance from the 
surface of neutron star because of the (1 — cos6 P7 ) _1 fac- 
tor (Link & Burgio 2005). Requiring conversion to take take 
place in the range R ^ r ^ 1.2R (at r — 1.2R the value of 
(1 — cosOp-,)^ 1 averaged over surface becomes double), P c has 
been found as ~ 0.021^ kEV (Link & Burgio 2005). From 
a follow up calculation by Link & Burgio (2006) it is found 
that conversion probability is slightly lower than that men- 
tioned above and can be parameterized as P c ~ Tq A kEV . 
Thus the total flux of neutrino/gamma ray generated in pul- 
sar from the decay of A + resonance is 

I = 2c£A pc f d (l - f d )n„P c (4) 

where £ is 4/3 and 2/3 for gamma rays and muon neu- 
trinos respectively. Denoting the duty cycle of the gamma 
ray/neutrino beam as ft (typically /& ~ 0.1 — 0.3), the phase 
averaged gamma ray /neutrino flux at Earth from a pulsar 
of distance d is given by 

0~ 2c£C»?/6/d(l -/<iH(5) 2 p c (5) 

where £ represents the effect due to neutrino oscillation 
(the decays of pions and their muon daughters results in 
initial flavor ratios 4> Ve : 4> v ^ : (f>v T of nearly 1:2:0 but at 
large distance from the source the flavor ratios is expected to 
become 1:1:1 due to maximal mixing of and v T ). C, = 1 
and 1/2 for gamma rays and muon neutrinos respectively. 

Average energy of the produced muon neutrinos would 
be e„ M ~ 50T-f keV TeV (Link & Burgio 2005, 2006) whereas 
that for gamma rays is expected to be E-y ~ lOOT^ 1 keV 
TeV. 



3 TEV GAMMA RAYS FROM FEW 

POTENTIAL PULSARS: COMPARISON 
WITH OBSERVATIONS 

Though there are about 1800 pulsars known through ra- 
dio detections, only few have been detected in the gamma- 
rays. From observations made with gamma ray telescopes 
on satellites so far only seven high confidence gamma-ray 
pulsars are known in the energy range up to few GeV 
(Thompson 2003). Besides, three other pulsars (PSR B1046- 
58, B0656+14, J0218+4232) are considered to be gamma- 
ray emitters with a lower confidence level (Thompson 2003) 
in the same energy region. Some basic properties of the high 
confidence young gamma-ray pulsars, namely distance (D), 
spin-down age, period (P), the calculated magnetic field 
strength at the neutron star surface, surface temperature 
and pulsar duty cycle (ft) are listed in the Table 1. 

None of the listed pulsars are, however, detected at TeV 
energies despite recent great improvement in the knowledge 
of Galactic gamma-ray sky above 100 GeV mostly by means 
of ground-based Imaging Atmospheric Cherenkov telescope 
systems such as the CANGAROO, HEGRA , High En- 
ergy Stereoscopic System (HESS) or the Major Atmospheric 
Gamma-ray Imaging Cherenkov Observatory (MAGIC). Up 
to now no evidence for pulsed emission from any other pulsar 



has been found from the observations (Chadwick et al. 2000; 
de Naurois et al. 2002; Lessard et al. 2000; Aharonian et al. 
2004; Aharonian et al. 2007), and only upper limits on the 
pulsed VHE gamma-ray flux are derived under various as- 
sumptions on the characteristics of the pulsed emission. For 
the pulsars listed in the table 1 the observed upper limits 
of integral fluxes are given in the table 2. The upper limits 
are obtained from the differential spectra assuming a power 
law differential spectrum with index —2.4. The upper limits 
of integral fluxes of gamma rays from the pulsars at around 
hundred TeV as obtained by the extensive air shower ex- 
periments such as the Tibet (Amenomori et al. 2008; Wang 
et al. 2008) are found less restrictive. 

Assuming pulsar accelerated ions are protons, numerical 
values of the integral TeV gamma ray fluxes are obtained for 
the pulsars listed in Table 2 from Eq.(5) and are also shown 
in the same Table (for the numerical estimation of flux we 
take Z = 1, and f d = 1/2 throughout the present work). 

It is clear from the Table 2 that the model with r\ — 1 is 
not consistent with the observed upper limits; the observed 
limits are substantially lower than the predicted fluxes. The 
observations made so far, however, do not rule out the model 
with r\ given by Eq.(3). 



4 TEV NEUTRINOS FROM PULSARS 

When r] is given by Eq.(3), neutrino fluxes from nearby 
young gamma ray pulsars would be much lower than es- 
timated by LB (Link & Burgio 2006). A high energy muon 
neutrino is usually detected indirectly, through observation 
of the secondary high energy muon produced by the muon 
neutrino on interaction with the ice or rock in the vicin- 
ity of a neutrino telescope via charged current interactions. 
The track of the produced muon is usually reconstructed 
by detecting the Cerenkov light emitted by the muon as 
it propagates through the telescope. The probability of de- 
tection of muon neutrinos is the product of the interaction 
probability of neutrinos and the range of the muon and is 
p„^ M ~ 1.3 x 10~ 6 (e„ M /l TeV) (Gaisser et al. 1995). The 
expected event rates in a neutrino telescope due to potential 
young pulsars are given in Table 3. 

The event rates are clearly very low and thus possibility 
of observing pulsars by a kilometer scale neutrino detector 
does not look bright. 



5 GAMMA RAYS AND NEUTRINOS FROM 
PULSAR NEBULA 

A young neutron star is generally encircled by pulsar wind 
nebula. Positive ions, after gaining energy from polar gaps 
will move away from the pulsar practically along the open 
field lines and will finally inject into the nebula. It is very 
likely that these energetic ions would be trapped by the mag- 
netic field of the nebula for a long period and consequently 
they should produce an appreciable of high energy gamma 
rays/neutrinos by interacting with the matter of the nebula. 
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Source 




d 

kpc 


age 
year 


V 

ms 


Bl2 


To.lkeV 


fb 


Crab 




2 


10 3 


33 


3.8 


~ 1.7 


0.14 


Vela 




0.29 


io 4 - 2 


89 


3.4 


0.6 


0.04 


B1706 - 


44 


1.8 


10 4.3 


102 


3.1 


1 ? 


0.13 


B1509 - 


58 


4.4 


10 3.2 


151 


0.26 


1 ? 


0.26 


J0205 + 


64 


3.2 


10 2.9 


65 


3.8 


0.04 


0.9 



Table 1. The characteristics of high confidence young low energy gamma ray pulsars 



Source 


Expected 


integral flux 


Observed upper limit of integral flux 




n = l 

10- 15 cm- 2 s- 1 


r\ as given by Eq. (3) 
lO-^cm- 2 ^- 1 


10- 15 cm,- 2 s- 1 


Crab 


1012 


1.65 


8 (56) 


Vela 


208 


0.123 


10 (20) 


B1509 - 58 


67 


.0034 


10 (20) 


B1706 - 44 


71 


.0025 


10 (20) 



Table 2. Comparison of predicted versus observed integral TeV gamma ray fluxes from some nearby young gamma ray pulsars. The 
numbers in par entheses are the energy thresholds in TeV for which upper limits are determined. The observe d upper limits fo r Crab 
pulsar is due to lAharonian et al.l l l2006al) whereas for the rest of the pulsars the observed upper limits are due to Aharonian ct al. | (120071 1 



Source 




Expected event rates 






km~ 2 yr~ 1 


Crab 




0.008 


Vela 




.0006 


B1509 - 


58 


.00016 


B1706 - 


44 


.00012 



Table 3. Expected event rates in a neutrino telescope due to some nearby young gamma ray pulsars. 
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5.1 Magnetic trapping of pulsar accelerated PeV 
ions in nebulae 

Conservation of magnetic flux across the light cylinder en- 
tails that outside the light cylinder B ~ r" 1 whereas far 
from the light cylinder radial component of magnetic field 
varies as B r ~ r~ 2 . Thus (far) outside the light cylinder the 
azimuthal component of the magnetic field dominates over 
the radial field. Therefore, accelerated protons while mov- 
ing away from the pulsar have to cross the field lines (for 
instance magnetic field lines at wind shock). The Larmor 
radius of particles (even for proton) of energy about 1 PeV 
is expected to be smaller than the radius of nebula during 
most of the time of the evolution of nebula (Bednarek & 
Prothcroe 2002; Bhadra 2006). Thus it is very likely that 
energetic particles of PeV energies would be trapped by the 
magnetic field of the nebula. The energetic particles propa- 
gate diffusively in the envelope and they escape from nebula 
when the mean radial distance traveled by the particles be- 
comes comparable with the radius of nebula at the time of 
escaping. This time is somewhat uncertain due to uncer- 
tainty of the value of diffusion coefficient but is estimated 
as at least few thousand years (Bednarek & Protheroe 2002; 
Bhadra 2006). 

5.2 Gamma rays and neutrinos from nebulae of 
young pulsars 

As pointed out in the preceding section the pulsar injected 
ions of PeV energies should be trapped by the magnetic 
field of the nebula for a long period and consequently there 
would be an accumulation of energetic ions in nebula. These 
energetic ions will interact with the matter of the nebula. 
The rate of interactions (£) would be nca p A, where n is the 
number density of protons in nebula and o p a is interaction 
cross-section. In each such interaction charged and neutral 
pions will be produced copiously. Subsequently the decays 
of neutral pions will produce gamma rays whereas charged 
pions and their muon daughters will give rise to neutrinos. 

If m is the mean multiplicity of charged particles in 
proton-ion interaction, then the flux of gamma rays at a 
distance d from the source roughly would be 

7 «2c / Sr 7 /„(l- /dR(f) irnt (6) 

where f3 represents the fraction of pulsar accelerated 
protons trapped in the nebula and t is the age of the pul- 
sar. Note that there should not be any reduction of flux 
due to pulsar duty cycle in the case of emission to nebula. 
Though n is taken as constant but actually at the early 
stages of pulsar n should be larger owing to the smaller 
pulsar period. So the above expression gives only a lower 
limit of flux. Typical energy of these resultant gamma rays 
would be ~ 10 3 /(6*m) TeV where for (laboratory) collision 
energy of 1 Pev m is about 32 (Alner et al. 1987). 

Numerical values of the integral TeV gamma ray fluxes 
from two nearby nebulae, Crab nebula and Vela nebula have 
been estimated for perfect trapping of pulsar accelerated 
protons in nebulae and are shown in Table 4. The observed 
integral gamma ray fluxes above 1 TeV from Crab nebula 
(Aharonian et al. 2006a) and Vela nebual (Aharonian et al. 
2006b) are also given there for comparison. 



The neutrino fluxes from the nebulae would be of nearly 
the same to those of gamma rays. Incorporating the neutrino 
oscillation effect the expected event rates in a neutrino tele- 
scope due to TeV muon neutrinos from nebulae of Crab and 
Vela are 0.2 km~ 2 yr~ 1 and 0.1 km~ 2 yr~ 1 respectively. Note 
that the event rates obtained here are rough numerical val- 
ues. The flux will be higher if the accelerated ion is heavier 
than proton. 



6 CONCLUSION 

To summarize, the present work suggests that pulsars are 
unlikely to be strong sources of TeV neutrinos. The non 
detection of any statistically significant excess from the di- 
rection of any pulsar by the AMANDA-II telescope (Acker- 
mann et al. 2008; Ackermann et al. 2005; Ahrens et al. 2004) 
is thus as per expectations. 

If protons are accelerated to PeV energies by the pul- 
sar, then pulsar nebulae are more probable sites of energetic 
neutrinos provided energetic particles of PeV energies are 
efficiently trapped by the magnetic field of the nebulae. But 
even for pulsar nebulae the expected event rates are small 
and the detection probability of pulsar nebulae by the up- 
coming neutrino telescopes, such as the IceCube (Halzen 
2006), is very low. 
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PWNe 


n 


Expected flux 


Observed flux 




cm~ 3 


xlQ- 12 cm- 2 s- 1 


10- 12 cm- 2 s- 1 


Crab nebula 


150 


0.6 


22.6 


Vela nebula 


1 


0.4 


12.8 
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